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In this paper we have analysed in detail two different purification protocols. The first one, 
proposed by Sudarshan, is based on the preservation of probabilities. We have constructed a second 
protocol here based on optimization of fidelities. We have considered both complete and partial 
measurements and have established bounds and inequalities for various fidelities. For every type of 
measurement, we have analysed post-measurement states based on the Maximum Entropy principle 
as well as what we have proposed as unbiased states. We show that our purification protocol always 
leads to better state reconstruction. These schemes can be thought of as operations in the sense of 
Kraus and we have explicitly constructed the Kraus operators for these. We have also shown that 
the entropy either increases or remains the same depending on the choice of the purification basis. 



I. INTRODUCTION 

In the density matrix formulation 1] the process of 
going from a mixed state to a pure state is called Pu- 
rification. There is a vast literature on this topic @. 
In the current literature 'Purification' is mostly under- 
stood to be the process of associating a suitable pure 
state of a larger system whose reduced density matrix 
is the mixed state one started with. This necessarily in- 
volves entanglement. We interpret purification in a larger 
sense to mean any protocol that produces a pure state 
from a mixed state. In this paper we consider purification 
schemes without entanglement. 

We specifically consider two schemes for purification 
the first of which was introduced by Sudarshan ^2] based 
on preservation of probabilities (scheme A), and a second 
proposed by us here based on optimal fidelities (scheme 
B). We apply both these schemes to the problem of recon- 
structing the pure state before a measurement from the 
post-measurement mixed state. We do so for both par- 
tial as well as complete measurements. In all such cases 
there is still a lot of freedom in constructing the post- 
measurement state itself. We have focussed on two spe- 
cific choices: The so-called maximum entropy state p ma x 
as well as an equal mixture p un b of the subensembles 
resulting from each measurement. We call the latter the 
'unbiased state' Q|. Operationally it is straightforward 
to realise as one has to simply divide the initial ensem- 
ble of pre-measurement state into equal parts for each 
measurement. These are discussed in Sec. IV. A. 

We apply our considerations to qubits only. We sep- 
arately consider the cases of a) complete measurement 
involving measurements of S x ,S y ,S Zl b) partial mea- 
surement involving measurements of S y , S z and finally c) 
which involves measuring S z only. For each of these cases 
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we apply the two purification schemes mentioned above 
to both p m ax , Punb ■ We show that in all cases scheme-B 
produces pure states with better overlap with the un- 
known initial pure state. 

We then show that there exist well defined operations 
in the sense of Kraus 0,0 which realise all these purifica- 
tion schemes. We explicitly construct the relevant Kraus 
operators. Though an 'environment' system is needed for 
realising these operations, our schemes do not require any 
entanglement between the qubit and the environment. 

We end our paper with an explicit proof that the total 
entropy of the qubit and the environment either remains 
the same or increases, in accordance with the second law. 
Though this is expected, we thought it desirable to ex- 
plicitly prove it. 



II. SOME PURIFICATION SCHEMES 

A. Purification Protocol - A 

Consider some density matrix which can be regarded 
as the mixture of two orthogonal states 



P = PiPi + P2P2 



(1) 



where, p{ = p\,p 2 = P2,tr(pip 2 ) = 0,trpi = trp 2 = 1. 

The Purification Protocol discussed here is based on 
the principle of preservation of probabilities. In 2] it was 
taken to mean that the overlap of the purified state p^ 
with p lj2 is pi j2 . 

This Protocol leads to the family of pure states: 



P 



(A) 



PlPl +P2P2 + \JP\P2 



/OlB> 2 + P2^P\ 



(2) 
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where II is a projection which is not orthogonal to either 
pi or p 2 . If pi = \0)(0\,p 2 = an d II is of the form 

IxKxl with x = p\0)+ v \\ 0»j v ^ since |x) is not 
orthogonal to either |0) or |1)) then the purified state is 
given by (4> is the phase of pv*) 



P2C 



1 



(3) 
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The reason that only this phase appears in the purified 
p is that preservation of probabilities leaves only a phase 
left unspecified in a pure state. Different choices of <f> 
lead to different purified states. There is no principle 
that selects a particular value of <fi. This protocol can 
only be implemented probabilistically. 

B. Purification Protocol - B 

Consider the state 

which we wish to 'purify'. The principle we adopt to fix 
the purified state is: the purified state must have maxi- 
mal overlap with the mixed state we started with. (The 
overlap between states p\ and pi is T — tr{p\.p2)). This 
amounts to demanding that the purified state be as close 
as possible to the mixed state. In contrast to A this is a 
deterministic protocol. 

It is elementary to show that the pure state p^ with 
maximum overlap with p is the eigenstate of p corre- 
sponding to the maximum eigenvalue . To see this let 
A+ be the maximum eigenvalue(considered degenerate 
for simplicity) and let the others be \,i = 1, N — I 
where N is the dimensionality of the Hilbert space, and 
let | A) be the corresponding eigenstates. Denoting the 
pure state by \ip) = a + |A + ) + J2i a i\K) we have 

i w = A + + ^H 2 (A J -A+) (5) 

i 

Since A.; — A + is negative it is obvious that the overlap 
is maximised when cti = i.e = |A+) and that the 
maximum value of the overlap equals A+. It is easy to 
show that 

A+ = a + A/2; A = I - 2a + y/4\p\ 2 + (I - 2a) 2 (6) 
and 

n 2 ( pA/2 \ (7) 

A(A + 2a-l) {p*A/2 A 2 /A J [ '> 

III. KRAUS FORMALISM AND 
PURIFICATION. 

Kraus et al 0, Q have given a formalism to study all 
possible changes of quantum states through the so-called 
operations. An Operation O is defined as follows: Con- 
sider a quantum system in the state p sys with a Hilbert 
Space TL which is coupled to another quantum system, 



often called the environment, in the state pe and which 
has a state space He- The system and the environment 
interact through a Unitary Evolution U which acts on 
the total Hilbert Space Tt ®He- Some property of the 
environment is selectively measured (what this means is 
that one selects a particular measurement outcome) , rep- 
resented by a projection operator Qe so that the com- 
bined state becomes: 

p= {l®Q E )U{p sys ® peW^I^Qe) (8) 

The system is then described by the reduced density ma- 
trix, p S y S — TrE p~, where the trace is taken over all 
possible states of the environment. The resulting state 
change O : p sys — > Psys is called an Operation. This op- 
eration can be represented in terms of Kraus operators 
Ak acting on the state space of the system such that 

Psys — ^ ^ AkPsysA^ (9) 
keK 

The operators Ak arc defined by 

(io) 

where /, g arbitrary vectors in the state space of the 
system, {fP \keK} are an orthonormal basis of Qe'H.e 
extended to "He and g E is the pure state in which the 
environment can be assumed to have started in. As the 
measurement Qe is selective, the Ak operators satisfy the 
trace non-increasing condition ^2 kfK AkA\ < I where K 
is some indexing set. 

Operations can connect any given pair of density ma- 
trices {pi,p2}, and in particular, a mixed state and a 
pure state. The entropy aspects of this are discussed in 
section V. 



A. Kraus Operators for Qubits 

We shall relax the condition of selectivity in measure- 
ments and consider all possible outcomes for measure- 
ments. We shall also restrict ourselves to qubits only. 
Then one needs two Kraus operators for a general oper- 
ation. 

Consider some orthonormal basis (called purification 
basis) |0), |1); and any pair of operators 

A = \yj)(0\ -Ai = 1^><1] (11) 

satisfying A^Aq + A\Ai = I. For any arbitrary density 
matrix pi„, these operators produce 

Pout = AopAl + A lP A{ = \ip)(tp\ (12) 

Clearly, p out is a pure state and it is independent of the 
initial state pi n . By ean (|f Ofl the Kraus operators here 
are of the form, 
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where the initial state of the environment is the pure 
state \0e)- It is straight forward to check that the unitary 
operator U that generates the Kraus operators of eanl jll[l 
is (with the notation \ip,x) — \tp)\x)E, (€\i/ J ) = 0): 

1 

U = 52(\1>,8){s,0\ + \t,s){s,l\) (14) 

s=0 

Kraus Operators relevant for the purification schemes of 
Sec. II can easily be worked out. 

IV. STATE RECONSTRUCTION FOR QUBITS 
AND EFFICIENCY OF PURIFICATIONS. 

Now we apply these considerations to address the is- 
sue of optimal reconstruction of unknown quantum states 
through both complete as well as partial measurements. 
The question we answer is: 'which of these schemes pro- 
duces a pure state that is closest to the pre-measurement 
pure state?' 

A. Post measurement state 

Clearly there are many ways of performing both com- 
plete as well as partial measurements leading to differ- 
ent post-measurement (mixed) states. Here we consider 
two particular ways of obtaining the post-measurement 
state: i) first is based on the well known maximum en- 
tropy principle and we shall denote the correspond- 
ing post-measurement state by p max - Such states have 
been studied extensively in [3j; ii) second is a proposal 
we are making in this paper which we call the unbiased 
mixture Q wherein the mixed states arising from differ- 
ent measurements are equally mixed to give the post- 
measurement state Punb- Since the first is discussed at 
length in Q we briefly describe the second. 

Consider an ensemble containing N identical copies 
of some unknown quantum state described by the pure 
density matrix /9j n j. One then subdivides this ensem- 
ble into k identical subensembles, makes measurement of 
S z (equivalently S y or S x ) if k = 1, makes measurements 
of S z and S y (equivalently S y and S x etc..) if k = 2, and 
finally measures all of S x , S y , S z if k = 3. One then puts 
together the resulting mixed ensembles in each case to 
form the post-measurement mixed state p un b- 

Now we construct these post-measurement states and 
apply the two purification schemes to each case. 

We use the eigenstates of a z , |±) z as the basis for the 
Hilbert space. The eigenstates of o~ x , o~ y are given by 
\±) x = 4-CI+), ± |-),); \±) y = M\+) z ± »!->,). 



B. Complete measurement: S x ,S y ,Sz measured. 

Let pi be the probability for the outcome |+) z , P2 for 
and P3 for |+)a;. The post-measurement density 



matrix for the S z measurement is: 

Pl=Pl\+)zz{+\ + {l-pi)\-)zz{-\ (15) 

and similarly for the other two subensembles. Now one 
takes an equal weightage of the three post-measurement 
density matrices to give p un b,3 = (l/ 3 )(/°i + P2 + Pa)- 
With the notation Aj = 2pj — l;j = 1, 2, 3 we can write, 

If Ai+2, A 3 -iA 2 \ f . 

P^ = e[A 3 +iA 2 3-Ai J (16) 

where '3' denotes the number of orthogonal spin compo- 
nents measured. Knowing pi,p 2 ,p 3 amounts to a com- 
plete measurement and the initial density matrix has to 
be 

If 1 + Ai A 3 -zA 2 \ 
Pini - 2 y A3 +iA 2 l-A x ) [U) 

Clearly, the relation between p U nb,3 and pm is, 

Punb,3 

= (1/3){I + Pim ) (18) 

In an earlier publication Q we have established the re- 
sult analogous to ean l|18|) for arbitrary systems with finite 
dimensional TL . The Fidelity of p un b,3 with the initial 
state is: 

F(Pini, Punb,a) = tr(p lm p unbi3 ) = 2/3 (19) 

independent of pi n i. Since, is pure, its eigen- values 
are 0, 1. Hence, the eigenvalues of p U nb,3 are, from 118fl . 
1/3 and 2/3. Therefore, 

Punb,s = l\l){l\ + l\s){s\ (20) 

Substituting in can ll8|) and using the completeness rela- 
tion I = \l)(l\ + \ s)(s\, one finds 

Pini = \l)(l\ (21) 

The eigenvector of p U nb,3 corresponding to the largest 
eigenvalue is the initial state. The Purification of l|2(J|) 
by Protocol - A is (see eqnlJSJ) 

Piil,3 - (2/3)|f><l| + (l/3)|*><*| 

+ V2/3(e^|0(s|+e-^| S )( /|) (22) 

This is clearly not the initial state. In fact, the fidelity 
of this state with the initial pure state is still only, 

npini,p[%)=m (23) 

On the other hand, the purification of p un b.3 by 
protocol- B is simply to select the state with the largest 
eigenvalue i.e. 

^,3 = |i><l| 



4 



which is indeed the initial pure state by virtue of 
can H21|) . The initial state is the closest pure state to the 
mixed state p un b,3 Protocol - B has purified the post- 
measurement state p un b,3 to maximum fidelity(in this 
case 100%) with the initial state. 

State reconstruction by the MaxEnt Principle 0, Q 
gives: 



Pmax.3 = + ^i ^ + A 2 a v + A 3 a x ) 



Pi 



(24) 



Clearly, in the case of a complete measurement state re- 
construction by the MaxEnt Principle and purification 
by protocol-B both yield the initial state. 

It should also be emphasised that even if one did not 
know that pt were the probabilities of outcome from com- 
plete measurements, but only knew the form of p un b,3, 
B would nevertheless reconstruct the initial state. The 
power of B is clearly evident in the case of partial mea- 
surements, discussed next . 



C. Partial Measurements: S y ,S z measured. 

In the case of partial measurements the initial state can 
never be unambiguously reconstructed. In this section we 
establish the following two results: (i) the purified state 
under protocol-B always has a greater fidelity with the 
initial state than does the post-measurement state, (ii) 
the fidelity of the purified state under protocol-B with 
the initial state is always greater than that of the puri- 
fied state under protocol-A (in an phase-average sense as 
protocol-A does not favour any single pure state) except 
in some singular cases where the fidelities are the same. 
Thus protocol-B is the better when trying to reconstruct 
the initial state from the post-measurement state. 

Suppose two measurements are made. Let p\ and p 2 
by the probabilities of outcome for \+) z and \ +) y respec- 
tively. Then post-measurement state (in the notation 
defined in Sec. IV. B) is, 



1 

Punb,2 — T 



Ai + 2 -iA 2 

iA 2 2 - Ai 



Let, the initial state be ipi. 



a 



01- 



(25) 



Then, 



(l+Ai)/2 = |a| , Im(aj3*) = A 2 /2. From these relations 
we can compute the fidelity: 



where A = \A\n. Withpi = + and p 2 = {\-^f), 
this is of the form (J). 

The initial state from l|17|) is, 

Pirn = \ + ± \ y/l - \A\ 2 a x (28) 

Denoting pi = |+')(+'| and p 2 = |— ')(— '\ the purification 
of p U nb,2 by protocol-A is, 



P { ul,2 = (5 + ^)\+')(+'\ + (5 " 



•1 , \A\ 



1 L4|. 



1 \A\ 2 

4~ 16 (l+')<-'K :/3 + |-')(+'| e -^)(29) 



Then 



tr{p un ba,Plnb,2) = 2 + ~§~ 

The eigenstates of p un b,2 are: 

|±') T = N±(iA 2 , A\ =F \A\) 



(30) 



(31) 



where N± 2 = 2{\A\ 2 t \A.A x |) are the normalization fac- 
tors. Using A+.A_ = 2\A.A 2 \ and {+'\cr x \-') = -ie iS ( 
5 is the phase of A 2 ) and (i'lc^i') = we can compute 
^(pinii Pmb 2 )- This is a function of the arbitrary phase 
P (from protocol-A). We therefore average the fidelity 
over the (3. This can be done in two ways: i) Average 
over all possible values of the f3 and ii) Average over only 
those values of the /3 (two such values) which gives opti- 
mal fidelity with the two possibilities for the initial state 
(eqnJ35|)). After some tedious work one gets 



J~av{,pinii P un b ,2) — 3~(Pini, Punb,2) 



(32) 



Using eqn(30) and \A\ < 1 one sees that this fidelity can 
never exceed 5/8 for any initial state. Now, it can be 
verified that purification by protocol-B (J7J gives, 



(B) 
Punb,2 



where 



P -ie tS ^fp{l-p) ■ 



Hpini,Punb,2) = -{A\+A 2 2 + 2) 



(26) 



where p ini = \ip ini )(ipini\- 

To purify the state by protocol - A, we can adopt the 
following procedure: From 125|) we find the eigenvalues of 
Punb.2 to be pi = (2+ \A\)/4 and p 2 = (2 - |A|)/4, where 
\A\ = (A 2 + Al) 1/2 . Using the notation A = (Ai,A 2 ) 
and a = (er z , a y ) we can write p un b,2 as a mixture of two 
orthogonal pure states: 



punb,2 = + 



n.cr. .1 

— )+P^2 



n.cr, 



(27) 



P= + 
Using the relations for |a| 2 , A\ and A 2 



HPini,P^b,2) = ^ + \A) 



(34) 



(35) 



This can indeed reach unity for a large class of initial 
states. From (l3*5|) . <|2lj|) it is clearly seen that, 

F{Pini, P u B n \. 2 ) ^ FiPini, Punb,%) (36) 
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State reconstruction by the MaxEnt Principle @] leads pj„, is of the form: 



to: 



_l(l + A x -iA 2 

Pmax.2 - I I- A X 



and we have the relation: 

Punb,2 = j (I + 2p TOa:Ci 2) 
Using this relation we get the fidelity: 

1 + \A\ 2 



3~ \Pini i Pmax ,2) — 



(37) 



(38) 



(39) 



Relation 138|) implies that the eigenstates of p m ax,2 are 
the same as that for p un b,2- Further, the eigenvalues can 
also be calculated directly from They are q\ — 

(1 + \A\)/2 and q 2 = (1 - |A|)/2. Purification of p m(KC ,2 
by protocol-A gives, 



(a) A + \A\ , , 1 + \A\ 

Pmax,2 = ( o )l+ )(+ I + ( o )|- )(- I 



4 



+')(-'|e i/3 +c.c.) 



(40) 



Again, using similar arguments as earlier in this section, 
we can calculate: 

•^aviPinij P max 2) — J~(Pinii Pmax,2 ) (41) 

Since the largest eigenvalue eigenstate of p ma x,2 is the 

same as t 
we have, 



same as that for p un b,2, we have p„f2c 2 = Pmb 2- Finally 



F{Pini, Pmax.2) ~ F{pini, P Mn b } 2) 
— J~ {Pini > Pmax ,2) = J~ av (pi n i , P max p) 
> F{pim,Punb,2) = Fav{pini, p[^{ n) ( 42 ) 



D. Partial Measurement: S z measured. 

When only one component of spin is measured, say, S z 
we have 

Punb.l =Pl\+)zz(+\ + (1 ~ Pl)\-)zz(-\ 

Again, the initial state is of the form: -0im = Q; |+)z + 
(3\— ) z . Now, we know only that \a\ 2 = pi. The Fidelity 
of the post-measurement state with the initial state is: 



F{pim, Punb,\) =Pi + (1 - Pi) 2 
Now, if we purify p un b,i by protocol- A, then 

Pi 



(A) 
PunbA 



VPi(l"Pi)e i0 



(43) 



Pi VpiC 1 ~Pi) e ~ 



(44) 



where p\ is known from measurement, but the phase 9 is 
unknown. The average Fidelity with equal weightage for 
all 9 is, 



Fav(pini, P„tLi) = Pi + C 1 - Pi) 2 



(45) 



Now, if pi > 1/2, then the purified state according to 
protocol-B is, 



(B) 
Punb 



(46) 



Then, 



Since, pi > 1/2, it can be verified that 

3~(piniTP un i l l) — 3~(PiniiPunb,l) — J~av(Pinii P un b l) 

' (48) 

It can be verified that even for pi < 1/2 protocol- B 
always leads to an improvement in fidelity. Thus, the 
fidelity offered by protocol-B is better than the average 
fidelity offered by protocol- A. 

In this case the state reconstruction by the MaxEnt 
Principle is: 



Pmax.l 



Pi 
1-pi 



Punb.l 



(49) 



Hence, the fidelity calculations remain the same as for 

Punb,l ■ 

V. ENTROPY OF PURIFICATION. 

As mentioned before, the purification process leads to 
interesting questions regarding entropy. Associated with 
every quantum state is the von Neumann entropy 



S v = —trplnp 



(50) 



In terms of the eigenvalues {A^} of the density matrix, 
one has 



S v 



^ A; lnA l 



(51) 



Hence the entropy S v is zero for pure states. On the 
other hand, for mixed states this entropy is always posi- 
tive. Thus purification reduces the entropy of the system. 
But the real issue is what purification, understood as an 
Operation on the composite of the system and the envi- 
ronment, does to the total entropy. 

To address this, let us first consider the action of the 
unitary operator U in eqn(|SJl which is given explicitly 
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in ean (|14fl . The initial state of the composite is p^ 
\0e)(0e\- K is easy to see that 



U(pW g, \0 E )(0 E \)Ui = \^)(i>\ <g> 



Pe 



(52) 



where p E is given by 



- ni°) 



Pe = P 



(0,0)|0 b )(0b|+p^(1,0)|1 s )(0b| 
-p^(0,l)\0 E ){l E \+p^(l,l)\l E ){l E \ (53) 



There are two remarkable properties of eans H52!53(l : i) 
the unitary transformation maintains the disentangle- 
ness of the system and environment and ii) the von Neu- 
mann entropy of the environment, which was zero to start 
with, has become equal to the system entropy we started 
with. So at this stage, not surprisingly, the total entropy 
has not changed(unitary transformations cannot change 
this!). This is an example of entropy swapping. 

Now we consider the effect of the projections Q E . Re- 
membering that for a total and not selective measure- 
ment we need two of them, the resulting environment 
density matrix is given by 



fin 

Pe 



p( o )(0,0)|0 £ )(0 £ |+p(°)(l,l)|l s )(l £ | (54) 



If the purifying basis |0), |1) is the same as the eigenbasis 
of , it is immediately obvious that the entropy of p E n 
is the same as that of p^ and the purification process 
does not change the total entropy. 

If that is not so, the entropy of p E n is still the entropy 
of a density matrix p^ which is diagonal and whose di- 
agonal elements are the same as that of p^ . It therefore 
follows that the determinant of p^ is greater than that 
ofp(°). Now we establish an important result: the von 
Neumann entropy of a two level system is an increasing 
function of the determinant of the density matrix. 

To prove this, note that the eigenvalues of the density 
matrix are given by A± = (1 ± ^/l — 4A p )/2 where A p 
is the determinant (0 < A p < 1/4). It is then straight 
forward to show that 



8Sy 

0A o 



1 



4A r , 



- 1 + >A^ >Q 



4A, 



(55) 



In this case the total entropy increases. 



VI. CONCLUSION 

In this paper we have addressed schemes for produc- 
ing pure states from mixed states. Though our schemes 
necessarily involve an environment, what is novel is that 
no entanglement between the system and environment is 
necessary. In this respect our purification schemes dif- 
fer from those that are mostly dealt with in the current 
literature. We give explicit expressions for the Kraus op- 
erators that are required. 

We have considered two purification schemes A and 
B as well as two post-measurement states, one of which 
is based on the maximum entropy principle. In all cases 
considered we find that scheme B purifies more efficiently 
than scheme A. The unbiased mixture state and the Max- 
Ent states perform equally well in all cases in producing 
a pure state as close as possible to the pre-measurement 
state. Operationally, unbiased mixture states are straight 
forward to realise compared to MaxEnt states. 

Finally we show that the total entropy either increases 
or remains the same. This is possible because we are able 
to achieve purifications without entanglement resulting 
in a case of entropy swapping. This along with a theorem 
we prove on entropy of two-state systems leads to the 
result. 
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